Invisibility and PT-symmetry 
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For a general complex scattering potential defined on a real line, we show that the equations 
governing invisibility of the potential are invariant under the combined action of parity and time- 
reversal (VT) transformation. We determine the "PT-symmetric as well as non-'PT-symmetric 
invisible configurations of an easily realizable exactly solvable model that consists of a two-layer 
planar slab consisting of optically active material. Our analysis shows that although PT-symmetry 
is neither necessary nor sufficient for the invisibility of a scattering potential, it plays an important 
role in the characterization of the invisible configurations. A byproduct of our investigation is the 
discovery of certain configurations of our model that are effectively reflectionless in a spectral range 
as wide as several hundred nanometers. 
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I. INTRODUCTION 

Consider a scattering problem for a real or complex potential v that is defined on the real line. The asymptotic 
solution of the time-independent Schrodinger equation has the form 
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where k is the wavenumber that takes real and positive values, and A± and B± are possibly fc-dcpcndcnt complex 
coefficients. We can obtain the information about the scattering properties of v using its transfer matrix M that is 
defined by the equation: 

r a i r a i 

(2) 

It is not difficult to show that M has a unit determinant This in turn implies that for every real or complex 
scattering potential the left and right transmission amplitudes coincide while this is not generally true for the left and 
right reflection amplitudes [![. Let T and R r ' 1 respectively denote the (complex) transmission and right/left reflection 
amplitudes, so that the transmission and right/left reflection coefficients are given by |T| 2 and |i? r,/ | 2 . We can express 
M in terms of T and R r ' 1 according to [|[ 
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Alternatively, we can describe the scattering properties of the system using its 5-matrix, 



S := 



T R r 
R l T 



(4) 



The potential is called reflectionless from the left (resp. right), if R l = and R r ^ (resp. R r — and R ^ 0.) As 
seen from (J3j) and ((4|), unidirectional reflectionlessness implies the non-diagonalizability of both M and S. Therefore, 
the parameters of the potential for which it becomes unidirectionally reflectionless correspond to exceptional points 
3 of M and S. 

The potential is called invisible from the left (resp. right), if it is reflectionless from left (resp. right) and in addition 
T = 1. Therefore, in light of we can express the condition for the invisibility of the potential from the left in 
terms of the entries of the transfer matrix according to 



Mm ^ 0, M 2 i = 0, Mn = M 22 = 1. 
Similarly, for invisibility from the right we have 

Mia = 0, Mn ^ 0, M u = M 22 = 1. 



(5) 
(6) 
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Because det(M) = 1, Equations ((5j) are not independent. The same holds for (j6)). We can reduce ([5]) and ([6]) to the 
following sets of equivalent independent equations respectively. 

M 12 + 0, M n = A/ 22 = 1. (7) 
Mai ^ 0, M n = M 22 = 1. (8) 

Recently Lin et al Q have reported the emergence of unidirectional invisibility for complex PT-symmetric locally 
periodic potentials of the form 

f a a + ae 2tl3z for \z\ < 
v(z) :=l 2 (9) 

\ for \z\ > f , 

where oto,a,/3,L are real parameters, and /? and L are positive. The spectral properties of the periodic potentials: 
v ( z ) = J2^Lo a n^ nZ , have been studied by Gasymov more than three decades ago @. The special case v(z) = ae % ^ z 
has been reexamined more recently in See also H. T he publication of 0] has motivated a more detailed study 

of the potential ©. In particular, the authors of [lCt use the analytic solution of the the wave equation 0, || 
to improve the approximate results of Q that rely on the rotating-wave approximation. For an earlier discussion of 
reflectionlessness and 'PT-symmetry, see fl2l |. 

In the present paper, we provide a thorough assessment of the role of PT-symmetry [l3T ] in the phenomenon of 
invisibility of scattering potentials in one dimension, and explore the invisible configurations of an exactly solvable 
model that admits experimental realizations. 

II. VT -SYMMETRIC NATURE OF INVISIBILITY 

We begin our analysis by examining the combined effect of parity and time reversal (VT) transformation on the 
transfer matrix M. It is easy to show that under this transformation the asymptotic solutions of the wave equation 
(fTJ) transform according to 

tp(z) ^ (VTtp)(z) := ip(-z)* -> A* ± e lkz + B* ± e- lkz for z -> Too, (10) 

where we suppose that k is real and positive. The transformation rule (|10[) implies that the transfer matrix of the 
PT-transformed system, that we denote by M^ 7 "*, satisfies 



A*_ 
B*_ 



A* 
B* 



(11) 



Combining Q and (jlip . we find the following transformation rule for the transfer matrix. 

M ^> M {vr) := A/T 1 *. (12) 

Because dctM = 1, this means that 

M U ^>M* 2 , M 12 ^-M* 2 , Afai^-Afa*!, M 22 ^ M n . (13) 

In particular, as noted in [3], the transfer matrix of a PT-symmetric potential satisfies M" 1 = M*. 

A straightforward consequence of (|13[) is that in general, regardless of whether the potential is PT-symmetric or 
not, the equations governing its invisibility (both from left and right), namely (0 and ©, are invariant under the 
PT-transformation |20j . We can restate this fact as follows. 

Invisibility Theorem: Consider a general real or complex scattering potential v that is defined on R. Let 
yCPi) t, e tfoe -pf -transform of v that is given by v^ p ' r ^(z) := v(—z)*, and k± be a positive real number. Then 
the following equivalent statements hold. 

i) v is invisible from the left (or right) for k — fc* if and only if so is v^^. 

ii) v is invisible from the left (resp. right) for k = fc* if and only if v* is invisible from the right (resp. left) 
for k = k*. 
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FIG. 1: (Color online) Schematic view of the cross-section of a two-layer infinite planar slab of gain material of thickness L 
that is aligned in the x-y plane. The arrows marked by R l and R r represent the reflected waves associated with incident waves 
from the left and right, respectively. 



For non-PT-symmctric potentials this theorem implies a pairing of the left/right-invisible configurations that are 
related by the P7~-transformation. We also have a pairing of all complex (non-real) potentials that display unidirec- 
tional invisibility at a given wavenumbcr where each member of the pair is the time-reversal (complex-conjugate) of 
the other. Obviously, parity reflection transformation changes a left- (resp. right-) invisible potential to a right- (resp. 
left-) invisible potential. In particular, even (P-symmetric) potentials do not support unidirectional invisibility. In 
view of the above theorem, the same holds for real (T-symmetric) potentials. Unlike P-symmctric and 7~-symmctric 
potentials, PT-symmetric potentials can display unidirectional invisibility. 

We can use the argument leading to the proof of the Invisibility Theorem to show that its statement also holds, if 
we replace the term "invisible" with "rcflectionless." In particular, we have PT-dual and T-dual pairings of complex 
reflcctionless potentials. 

III. TWO-LAYER INFINITE SLAB MODEL 

Consider an infinite planar two-layer slab of optically active material as depicted in Figure [TJ Suppose that the 
layers have equal thickness, L/2, and different gain/loss properties. Then the wave equation for a linearly polarized 
time-harmonic electromagnetic wave that propagates along the normal direction to the slab can be reduced to the 
Hclmholtz equation [15j . 

■$>"(z) + k 2 n{z) 2 ^(z) = {) 1 (14) 

where 

{tli for -f < z < 0, 
n 2 for < z < f, (15) 
1 for | Z |>§, 

and tti and n 2 are complex refractive indices of the layers. In Ref. [l5|, we have used this model to examine the 
importance of PT-symmetry for generating self-dual spectral singularities, i.e., spectral singularities Q that coincide 
with their time-reversal dual. The concept of a self-dual spectral singularity provides the mathematical basis for an 
optical device that functions as a coherent perfect absorber (CPA) [lj,[l6[ whenever it is subject to incidence coherent 
waves with identical amplitude and phase from both sides and operators as a laser otherwise. It is therefore called a 
CPA-laser [l7|- In we establish the existence of non-PT-symmetric CPA-lasers. Here we use the optical system 
defined by (|15|) to examine the role of PT-symmetry in achieving unidirectional invisibility. 
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We can calculate the transfer matrix M of the above model using the method outlined in [18J and employed in [15j . 
Doing the necessary calculations this yields the following expressions for the entries of M. 

n_n_sina_)}e-^, (16) 

n_n+sina+)j, (17) 

n_n+sina+)j, (18) 

(19) 



M u 


= « 


-n 2 _ 


) -1 |tt+ cos a + - n 2 _ 


cos a_ + 


i(n+n + sin a + 


M 12 


- (4 


-n 2 _ 


) _1 |n_n + (cos a + — 


cosa_) - 


f i(n + n_ sina_ 
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)~ 1 {rt 2 h cosa + - n 2 _ 


cosa_ - 
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n±:=ni±n 2 , a± := , K:—Lk, ft± := mtl2 ± 1 = ^(tt+ — tl?_) ± 1. (20) 



where 



As a nontrivial check on the validity of equations (fl6|) - (|T9| . we have shown by direct calculations that indeed 
detM = 1. 

In view of ([17]) - (|19|) , the condition ((5j) for the invisibility of the system from the left takes the form: 

n_n+(cos a+ — cos a-) + i(n+n_ sin a- + n_n + sina+) ^ 0, (21) 

n_n+(cos a + — cosa_) — z(n + n_ sin a_ + n_n + sino + ) = 0, (22) 

n+ cos a+ — n 2 . cos a_ + i(h + n + sin a+ + n_n_ sina_) = (n^ — n^_)e I,R , (23) 

x\ 2 + cos a + — n 2 _ cos a_ — i(n + n + sin a + + n_n„ sino_) = {n 2 + — X?_)e~ lSi . (24) 

Note that for ni = n 2 . n_ = and we can use Equations (|22|) and (|24|) to show that in this case tli = 1. This 
corresponds to an empty slab. Therefore, without loss of generality, we take n 2 ^ tii. Moreover, by adding and 
subtracting Equations ([24]) and ([23]) . we can reduce them to 

tt^_ cos o + — n 2 , cos o_ = (n + — n_) cos .ft, (25) 
n+n+ sina+ + fx n sin a_ = (n+ — n 2 .) sin.ft. (26) 

Similarly using (f2Tj) and (|22|). we have 

n_n + (cos a + — cos o_) ^ 0, (27) 
ti + n_sina_ + n_n + sina + ^ 0. (28) 

For non-exodic material, Re(tti i 2) > 1 and n+ ^ 0. We also have n_ ^ 0. Therefore, (f!?7| holds if and only if 

cos o+ ^ cos a_ . (29) 

We arrive at the same conclusion if we demand that the potential be invisible from the right. Therefore, a necessary 
and sufficient condition for the unidirectional invisibility of our system is that we satisfy (|2"5j) . and ([2T)|) . To 

determine the direction of invisibility we need to check which of M± 2 and M 2 \ vanishes. 



IV. BIDIRECTIONAL INVISIBILITY 



Suppose that ni ^ n 2 and the system is invisible from both directions, i.e., it satisfies ([25]). (|26|) . and cos o+ = cos a_. 
Substituting this relation in (f2"2")l and ([231) an( i noting that - n 2 ^ gives 

coso± = cos .ft, (30) 

rL|_n_sina_ +n_n+sino+ = 0, (31) 

Equation ([3"0[) means that 

a+ = 27TTO+ ± .ft, a+ = 27rm_ ± .ft, (32) 

where m± are integers. In view of ([20)1 . this implies that n±, h±, tli, and n 2 must take real values. Another consequence 
of (001) is 

sin a + = ± sin a_ . (33) 
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Inserting this relation in (|22j) and (|24[) and using the fact that n±, ri±, and a± are real, we find 

(n+n_ ± fi_n+) sin a + = 0, (n+n+ ± n_n— ) sin a + = 0. 

Adding these two equations gives 2rii(n + ± ri_) sin a + = 0. Because tii(n+ ± n_) 7^ and ([33]) holds, this implies that 
sin o± = 0. Equivalently, there are integers m'± such that 

o± = 7rm' ± , (34) 

Next, we combine this relation and (|32[) to obtain 

H = nm, (35) 

where m is an integer. In light of (f3"4"| and (|55|) . we also find that m' + ± m'_ are even, and 

2mi 2m 2 2L 

ni = , n 2 = , A= — , (36) 

m mm 

where m\ := (m+ + m'_)/2, 7712 := (m^_ — m'_)/2, and A := 2n/k = 2irL/& is the wavelength. 

In summary, our device is invisible from both directions provided that the refractive indices of both the layers are 
rational numbers, and the wave length A at which the device is invisible is such that L is a half integer multiple of A, 



2 lj| . Notice that even bidirectional invisibility does not require PT-symmetry. 



V. UNIDIRECTIONAL INVISIBILITY AND VT -SYMMETRY 

Equations (f25f and ([26| are complex equations involving two complex variables, n±, and a real variable A. Therefore, 
in principle, we can fix say the real part of n+ , which is bounded from below by 2, and try to solve these equations 
for the imaginary part of n+, the real and imaginary parts of n_, and A. For the physically relevant range of the 
values of these variables this turns out to be difficult to implement directly. This calls for a more systematic study 
of the structure of (|25p and (|26|) . An important clue is provided by PT-symmetry of the invisibility equations that 
we established for general scattering potentials in Section II. Here we offer a direct and explicit verification for the 
presence of this symmetry. 

V 

First we note that for our system the parity transformation means swapping the labels of tij, i.e., ni < — > 112, and 

T 

the time-reversal transformation corresponds to complex conjugation of the refractive indices, xij — > n*. Therefore, 

VT 

ni < — > rij, and consequently 

VT * VT * VT * VT * - VT ~* /0 „n 

n + — > n + , n_ — > — ti_, a + — > a + , a_ — > — a_, n± — > n ± . (37) 

Performing these transformation on Equations (|25[) and (|26p and noting that their left-hand sides are analytic functions 
of n±,a±, and n±, we see that 'PT-transformation maps these equations to their complex-conjugate. Therefore, the 
real and imaginary parts of the transformed equations are identical with those of (|25p and (|26p . This confirms our 
general result ensuing the PT-symmetry of the invisibility equations. The presence of this symmetry suggests that 
we first examine the PT-symmetric solutions of (f2~5j) and (|2"6")l . 



A. PT-Symmetric Invisible Configurations 

The optical system described by (fl~5j) is PT-symmctric provided that n| = ri2. In this case, n+ and n_ are 
respectively real and imaginary. The same is true about a±. In particular, if we denote by ry and k the real and 
imaginary parts of ni , we have 



2r?, 



2in, 



n±=rf + n 2 ±l = |ni| 2 ±l. 



(38) 



Using these relations and the identities cos(iz) = cosh z and sin(iz) = isinhz, we can express (|25[) and (|26p in the 
form 



rj 2 + k 2 
1 



cos (£77) 

1 



Tj 2 + K z 



rj 2 + k 2 



i] sin(^?7) — 11 — 



cosh(^K) = cos^, 
1 



I] 2 + K 2 



k sinh(^K) 



= sin^. 



(39) 
(40) 
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It is easy to see that (|29[) holds automatically. Notice that (|39|) and (|40)) are real equations involving three real 
variables, 77, k, and This is an enormous simplification over the general non-T'T-symmetric case that amounts to 
solving four real equations for five real unknowns. 

Next, we recall that for typical optically active material, |k| < 10 -3 and 1 < rj < 5. Therefore K 2 /rj 2 < 10~ 6 , and 
we can safely approximate rj 2 + n 2 by r] 2 . This simplifies (f39|) and (|4"0)) . and we find 



, „ . k 2 cosh(^ft) „ , , , 

cos(&77) + ^ -«cos£, (41) 

77 

(1 + rf 2 ) r) sm{Sbj) - (l - 77~ 2 ) « sinh(^K) sinii. (42) 

Here and in what follows we use the symbol "rj" to mean that we ignore terms of order k 2 /rj 2 . 

Because & = Lk = 2irL/\, where A is the wavelength, & can take very large values. However, according to (|39l) . 
cosh(^K) < 2t] 2 /k 2 . This in turn implies 

*<tWSV (43) 



Therefore the smaller |k| (or the gain/loss coefficient) is, the larger the upper bound on & becomes. In particular, (|43|) 
does not impose any severe restriction on the wavelength at which the system is invisible from the left. For example, 
using the fact that usually 1 < rj < 5 and |k| < 10 -3 , we find from (|4*5)) that 

R < 18420. (44) 

This is equivalent to L < 2930A. In particular, for A > 100 nm, we have L < 293 /urn, which can be easily realized 
cxp er imcnt ally. 

Another simple consequence of (|4"Tj) is that r\ = 1 implies k = 0. Again this corresponds to an empty slab. Therefore 
unlike for the 'PT-symmetric periodic potential considered in Q , the unidirectional invisibility can be achieved in our 
system provided that we use an optically active medium with 77 > 1 . 

Next, we obtain an analytic solution of (|41|) and (|42|) . To do this we introduce 



A:=(^) 2 , «:=£±±, /3 := (45) 
rj z — 1 rj z — 1 

respectively solve (|4Tj) and (|42|) for cosh(^.K) and sinh(^K), and employ the identity cosh 2 (^K) — sinh 2 (^K) = 1. This 
gives a quadratic equation in X with a single positive solution. We can use this solution to express k in terms of 77 
and & as follows. 



/ 1 4 2 1 2 

± 77 y y — [a sin(^r/) — /3 sin + [ cos(.S?7) — cos j?] — - [a sin(^r/) — f3 sin . (46) 

Using this relation in (|41[) we find an equation for 77 and ^. Because (|41[) involves cosh with possibly large arguments, 
it is not suitable for numerical calculations. Therefore, we solve this equation for cosh(.&c) and employ the identity 
cosh - (z) = ln(z ± \J z 2 — 1) to express it as 



£|re|-lnj ^ ( cos^-cos(^) ± y [cos£- cos(.&7)] 2 - ^ ] ^ « 0. (47) 




Now, we substitute (|4*5|) in the left-hand side of (|4"T)) . The result is a function of 77 and & that we denote by F±(rj,&). 
In this way, we can express (|4?p in the form 

F±fe^«0. (48) 

We can choose some typical values for 77 and solve (|4"5)) numerically to find the possible values of Once these 
are determined we use (|46[) to obtain the corresponding values of k. In order for the latter to be experimentally 
meaningful, they should be sufficiently small. We expect this to be the case provided that the & values we find fulfil 
the condition (|44[) . 

To be specific, suppose that the left-hand (right-hand) layer of the slab is the one containing the gain (lossy) 
material, so that k < 0. In this case, we find that F_(r),&) 7^ 0, while for each 77 > the equation F+(r],R) = gives 
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a (countably) infinite set of values for .ft. Using these values in (|46|) we can determine the corresponding values of n. 
The & and K values we obtain in this way satisfy (|48p . but they may not satisfy (|4"Tj) and (j4"2"|) . This is because in 
deriving (j48|) we affected functions on our expressions that arc not one-to-one. Therefore, although every solution of 
(|4"Tj) and (|4"2")l is a solution of (|4l)|) and (|4"8|) , the opposite is generally not true. In our numerical calculation, we found 
solutions of (|46f and ([48]) and checked if they satisfy the exact equations (|39| and (|40|) . In this way we determined 
the acceptable approximate solutions of the problem that are valid whenever we can neglect k 2 /r/ 2 . We then obtained 
the exact numerical solutions of and (|4"U1) in the neighborhood of the approximate solutions. For example for 
rj = 1.4,2.0,3.4 with & ranging over (1995,2005) we respectively found two, three, and five different values of X/L 
(between 317.7 and 318.9) for which the system displayed unidirectional invisibility. Except for a single data point, 
the approximate and exact results agreed to eight significant figures. We also checked the values of M 12 and M 2 2 to 
determine the directional of invisibility. It turns out that for the case that k < (resp. k > 0) we have invisibility 
from the left (resp. right). 

Note also that because the equations we used to characterize invisibility of our system are even in k, changing the 
sign of the values of k that lead for invisibility from the left, we find the values that support invisibility from the 
right . This confirms the statement ii of the Invisibility Theorem. 

Figure [5] shows plots of \T\ — 1, \R l \, and the argument (phase angle) of T as a function of M. for a PT-symmetric 
sample with r\ = 3.4 and k = —.00342163. The fact that the graphs have a common intersection point located on 
the M.-axis is a clear demonstration of the invisibility of the device from the left for the value of & at this point, i.e., 
J? = 2000.147552. For these values of 77, k, and .ft, we have 

|i?'| 2 < 10" 10 , |T| 2 -1<10- 5 , |iT| 2 > 0.89. (49) 

Taking the thickness of each layer of the slab to be 150 /im, so that L = 300 /Ltm, and using .ft = 2000.147552, we 
find that the device displays invisibility from the left for A = 942.408269 nm. If we keep only three decimal places in 



0.3 

0.2 

0.1 

0.0 
-o.i 
-0.2 
-0.3 

2000.05 2000.10 2000.15 2000.20 2000.25 
% 

FIG. 2: (Color online) Graphs of |T| 2 — 1 (thick solid curve), argument (phase angle) of T (dashed curve), and \R l \ 2 (thin solid 
red curve) as functions of J? = Lk = 2nL/X for a 'PT-symmetric two-layer slab with m = n| = 3.4 — 0.003422i. The intersection 
points of these curves and the .S-axis corresponds to invisibility from the left for A = 2000.148. For a sample consisting of two 
150 fim thick layers, this gives the value 942.408 nm for the wavelength. 

the above figures for k, M. and A, the first two bounds given in (|39"]1 change to |i?'| 2 < 10~ 5 and \T\ 2 - 1 < 3 x 10~ 3 , 
respectively. 

B. Non-PT-Symmetric Invisible Configurations 

In the absence of "PT-symmetry, the equations governing unidirectional invisibility, (|25[) and (|26p. are complex 
equations involving two complex and one real unknowns, n+, n_, and ^, respectively. It is more convenient to 
parameterize them in terms of a+, a_, and ^. This leads to 

a 2 + (cos a + — cos M) — o_ (cos a_ — cos &) = 0, (50) 
(a 2 + -a 2 _- & 2 )a+ sina+ + {a\ - a 2 _ +£ 2 )a_ sincu - 2(a 2 + - a 2 _)ilsinJ? = 0. (51) 
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Equating the real and imaginary parts of the left-hand side of these equations to zero, we find a set of extremely 
complicated real equations for & and 

x± := Re(a±), y± := Im(a±). (52) 

A direct Mathematica-based numerical treatment of these equations proves to be intractable for the physically relevant 
ranges of the values of x± , and y± . A careful study of the reasons for the failure of such a treatment guided us 
to a particular change of variables that allowed for a reliable pcrturbativc method of solving these equations. The 
following is a brief outline of this method. 

1. Because \x±\ = | Re(^n±)|/2 and & are typically large numbers, we express them as 

x± = 2irm± + — — , & = 27rmo + - — - — , (53) 

27rm± 27rmo 

where m± and mo are integers, m + and mo are positive, and j± and 70 are real parameters that are of the 
same order of magnitude as y±. In particular, |7±|, I70I and \y±\ are at least two orders of magnitude smaller 
than m_|_, |m_|, and mo- In view of (|53p and the fact that 

a±=x±+iy±, (54) 

this justifies the following approximations 

sinh y± 

coso± w coshy± , (55) 

27rm± 

a± sin o± « 7± cosh y± — y± sinh y± + r Litim± sinh y± , (56) 
0j_ « 47rm±(7rm± + iy±), cos£~l, AsmM.~jo. (57) 

Note also that because a± = &n±/2 = M.(rii ± n_)/2 and Re(tii) > 1 and Rc(ri2) > 1, Equations (|53[) and 
imply that 

|m_| < m+, mo < m+. (58) 

2. Substituting ((55]) - (|57|) in (I5TJ1) and (|5Tj) and keeping only the dominant terms, we respectively find from the 
real part of (|50p and the imaginary part of ()51[) the following approximate equations. 

coshy_ rj /i 2 (cosh2/+ — 1) + 1, (59) 
sinh?/_ rj — /i^sinhy + , (60) 



where 



:= — . " := —2 — —2 — =t- ( 61 ) 



m_ — m_ — m 2 . 



Similarly implementing the approximations (|55p - (|57p in the imaginary part of (|50[) and the real part of (|5T 
and using (|59[) and (|60[) to simplify the resulting equations, we obtain 



where 



M7+ sin hy+ -7-sinhy_ w 2/i(y+ - juy_)(coshy+ - 1), (62) 
7+ cosh y+ + 7_ cosh y_ w [(1/ + v + )y + - v yJ) sinh y + 

+ Woy+ + (1 - "-)2/-] sinhy_ + {y+ - vJ)^ Q , (63) 



2ml 2m_?7ii 



6 ± 

■ 9 o 95 ^0 • 9 9 9 ■ 

m^_ — m_ — mg m^_ — mi. — mg 

3. Inserting (f5U)) and in the identity cosh 2 y_ — sinh 2 y_ = 1 and using cosh 2 y + — sinh 2 y + = 1 to simplify the 
result, we find a quadratic equation in coshy + with solutions 

coshy + w 1, (64) 

coshy+ rj = — = - 65) 

M — v 1 - S^r 
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In view of (f59|) . (|64|) implies y± « 0. This corresponds to the bidirectional invisibility that we examined in 
Section IV. (p5|) is acceptable provided that its right-hand side is greater than 1. This leads to the condition: 
v 2 



,2 < /i 2 which, in light of ([58)) and (|6"Tj) . is equivalent to 



mo < m + — |m_|. (66) 
Supposing that this condition is fulfilled, we can invert (|65p to obtain 

fu 2 + v 2 -2\ 

y + « iarccosh - — . (67) 

V /x 2 - v 1 ) 



Furthermore, using (|60[) . (|65|) . and the identity: sinhy + = ±i/l — cosh y + , we find [2 



y_ » Tarcsinh I * " 1 ■ ( 68 ) 

Finally, we substitute ([57]) and (f6"51) in (IB"2l and (j6"3")l to obtain a pair of linear equations for 7± and 70 . We can 
easily solve any two of these variables in terms of the third. 

As we see our method includes four free parameters, namely m±, mo, and any one of 7± and 70. We can easily 
relate m± to the real part of ni and ri2, that we respectively denote by r/i and 772. More specifically, 

m ± « fa± 2 ??2)mo . (69) 
The integer mo is related to R = Lfc = 27ri / A according to 

mo « I = |. (70) 

Therefore, for a given two-layer slab and any wavelength A we can use the real part of the refractive index to determine 
m± and mo. The choice of the last free parameter, i.e., one of 7± and 70 is arbitrary. We choose 70 to be a real 
number of order 1 so that the corresponding 7± are also of order 1, [23|. It turns out that the choice of the + or — 
sign in (j6"T|) corresponds to invisibility from the right and left, respectively. 

The following is the result of the application of our method for a sample with 771 ps 3.4, ijj w 1.4, I « 2000, and 
70 = -6. 

m_ = m = 318, m+ = 764, R = 1998.049925, (71) 

y_ = 2.392197, y+ = -1.180878, (72) 

7_ = -1.387129, 7+ = 1.062178, (73) 
m = 3.402510 + i(6.062508 x 10" 4 ), n 2 = 1.402514 - 788281 x 10~ 3 ). (74) 

For L = 300 /im, this sample displays invisibility from the left at A = 943.397644 nm. Inserting the above numerical 
values in the exact expressions for the reflection and transmission coefficients gives 

||T| 2 - 1| < 2.1 x 10~ 5 , |arg(T)| < 3.2 x 10~ 3 , |i?'| 2 < 2.8 x 10~ 6 , |i? r | 2 > 14.1. 

Therefore our device is invisible from the left. The large value of |i? r | 2 indicates that it reflects the incident waves 
of this wavelength from the right after a 14-fold amplification. Furthermore, choosing the opposite sign for y± as in 
(1721 gives the time-reversed configuration: 

m = 3.402510 - i(6.062508 x 10" 4 ), n 2 = 1.402514 + £(1.788281 x 10~ 3 ), (75) 

that is invisible from the right and has a large reflection coefficient from the left (that is also greater than 14.). This 
is a manifestation of the statement ii of the Invisibility Theorem. 

Figure [3] shows the plots of the |i?'| 2 , \T\ 2 — 1, and arg(T) as functions of A and log|i? ; | 2 as a function of A for 
the non-'PT-symmetric sample with specifications (|74|) . An interesting outcome of our findings is that the reflection 
coefficient from the left remains negligibly small for a very wide spectral range. Although the invisibility from the left 
is lost for wavelengths slightly different from A = 943.397644 nm, the system remains essentially rcflcctionlcss from 
the left (with |i?'| 2 < 10 -4 ) for wavelengths ranging from about 860 nm to 1060 nm. 
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FIG. 3: (Color online) On the left are graphs of |T| 2 — 1 (solid blue curve), arg(T) (dashed purple curve) , and |7?'| 2 (thick solid red 
curve) as functions of J? for a non-PT-symmetric two-layer slab with ni = 3.403+i(6.063x 10 -4 ) and ri2 = 1.403 — i(1.788x 10 -3 ). 
The intersection points of these curves and the .ft-axis corresponds to invisibility from the left for R = 1998.050. For L = 300 /im, 
this gives A = A* := 943.398 nm. On the right is the graph of log \R l \ 2 as a function of the wavelength for L = 300 /im. The 
dashed line represents A = A» at which the system is invisible from the left. 



Next, we give the results for a sample with a different value of r/2- We take rji sa 3.4, r\2 ~ 2.0, M. w 2000, and 
7o = —2. These lead to 

m = 318, m- = 223, m+ = 859, R = 1998.051927, 

y- = 1.704272, y+ = -0.492952, 7 _ = 0.612879, 7+ = -1.307241, 

tii = 3.402514 + i(6.062502 x 10" 4 ), n 2 = 1.999999 - 099683 x 10" 3 ), (76) 
||T| 2 - l| < 2.4 x 10" 6 , |arg(T)| < 4.0 x 10~ 4 , |i?'| 2 < 7.0 x 10~ 7 , |i? r | 2 > 0.9. 

For L = 300 /im this configuration is invisible from the left at A = 943.396699 nm. Again changing the sign of the 
imaginary parts of tii and n-2, we find the dual configuration that is invisible from the right. Furthermore, for the 
original system (|T6"|) . the reflection coefficient from the left remains negligibly small (less than 10 -4 ) for wavelengths 
between 780 nm and 1165 nm. Therefore, the device is essentially reflectionless in a spectral range that is 385 nm 
wide. 

Finally, we should point out that if we keep only three decimal places in the numerical values we use in our 
calculations, the upper bounds we obtain for ||T| 2 — l|, |arg(T)|, and |i?'| 2 arc increased by about two orders of 
magnitude. 



VI. CONCLUDING REMARKS 



The search for optical realizations of 'PT-symmetric potentials have unraveled unusual and interesting physical 
properties 0, Among these is the unidirectional invisibility achieved for the PT-symmetric locally periodic 

potential given by 0, 3- I n this paper, we showed that the phenomenon of the invisibility of a complex scattering 
potential is PT-symmetric in nature. By this we mean that irrespective of whether the potential possesses VT- 
symmetry or not the equations characterizing invisible potentials are VT- invariant. This does not means that non- 
■pT-symmetric configurations cannot display unidirectional invisibility. It only means that the PT-symmetric invisible 
configurations are quite special, for they possess the same symmetry as the equations. A simple consequence of the 
presence of this symmetry is that the non-PT-symmetric configurations come in pairs that are related by the PT- 
transformation, i.e., they are PT-dual of one another. Furthermore, if a one-dimensional scattering system is invisible 
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from the left for some value of the wavenumber k, then its time-reversal will be invisible from the right for the same 
value of k. These VT- and T-dual pairings of invisible potentials also apply to reflectionless potentials. 

As a concrete implementation of our general results we have examined a simple two-layer planar slab of optically 
active material. We have studied in great detail the structure of the PT-symmetric as well as non-PT-symmetric 
invisible configurations of this model using essentially analytic techniques. In particular, we have shown by explicit 
calculations how "PT-symmetry simplifies the characterization of the invisibility of our model. We have also confirmed 
the above-noted T-dual paring of invisible configurations by showing that the time-reversal of various left-invisible 
configurations arc right-invisible. Furthermore, we found invisible configurations that remain reflectionless within a 
very wide spectral range. We intend to explore this phenomenon as a separate research project. 
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